Oscillatory dynamics in nano-cavities with non-instantaneous Kerr response 
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We investigate the impact of a finite response time of Kerr nonlinearities over the onset of spon- 
taneous oscillations (self-pulsing) occurring in a nanocavity. The complete characterization of the 
underlying Hopf bifurcation in the full parameter space allows us to show the existence of a critical 
value of the response time and to envisage different regimes of competition with bistability. The 
transition from a stable oscillatory state to chaos is found to occur only in cavities which are detuned 
far off-resonance, which turns out to be mutually exclusive with the region where the cavity can 
operate as a bistable switch. 
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I. INTRODUCTION 

Self-pulsing (SP), the onset of spontaneous oscillations, 
is a universal feature of nonlinear structures with feed- 
back. As long as passive systems are concerned SP has 
been investigated theoretically in settings ranging from 
isolated ring cavities [THS] and parametric intracavity 
mixing [3H6] to Bragg gratings [TJ [8] or grating-assisted 
backward frequency conversion schemes [9j [10] , and it is 
still a subject of active research [TTHT5] . In particular, 
the dynamics of nonlinear passive cavities, whose study 
has been pioneered in the eighties [IH5] , is extremely rich 
encompassing stable as well as chaotic SP which can com- 
pete with bistabilities and transverse effects [14]. His- 
torically, SP and chaos (the optical equivalent of strong 
turbulence) have been first analyzed by means of delay- 
differential models accounting for the round-trip delay at 
each passage, which can be large in ring cavities. An 
instability named after Ikeda occurs in this framework 
when the relaxation time of the nonlinear response is 
much shorter than the transit time [1 and has been 
tested experimentally [15]. However SP and so-called 
weak turbulence occurs also in the opposite (say short 
cavity) limit, where the delay can be averaged out to end 
up with a differential mean- field model [2 . This regime 
becomes important nowadays where nano-cavities are 
employed for many modern photonic applications [16] , in- 
cluding bistability [17 , demonstrated in photonic crystal 
(PhC) membranes which offer great flexibility of design 
as well as high nonlinear performances in semiconduc- 
tors [TM25] . In these cavities, transverse effects are ab- 
sent and their dimensions are so small that the response 
time of the medium can be much larger than the light 
transit time in the cavity, yet being comparable with the 
photon lifetime which is strongly enhanced on account of 
the large quality factor Q. SP in such nano-cavities has 



been recently predicted, owing to the free-carrier disper- 
sion induced by two-photon absorption [13]. The main 
features of such mechanism is the existence of a critical 
value of the time relaxation constant r, as well as a wide 
region of the parameter space where stable (non-chaotic) 
SP can be potentially observed. In this paper we analyse 
the dynamics of a nano-resonator in the good-cavity limit 
when the underlying nonlinear mechanism is a Kerr-like 
nonlinearity with finite response time. Our analysis is 
based on the differential model proposed in Ref. [2]. In 
spite of the simplicity of such model, which makes it an 
ideal prototype for understanding the role of relaxation 
processes, a full characterization of SP and its competi- 
tion with bistability in the full parameter space was never 
reported (to the best of our knowledge) after Ref. [2 . We 
propose such analysis, adopting a different normalization 
with respect to that employed in Ref. [2 , better aimed 
at capturing the key role of the relaxation time. This is 
especially important nowadays in view of assessing how 
a given designed nano-cavity may be expected to behave 
by changing the characteristic relaxation time of the non- 
linearity as a consequence of choosing different materi- 
als and/or adopting techniques for fine tuning their re- 
sponse time. We propose an analytical characterization 
of the SP instability and its competition with bistability 
in the full parameter space, pointing out the qualitative 
similarity with the features observed for a nonlinearity 
dominated by free-carrier dispersion [15] . Nonetheless, 
we further investigate also the destabilization mechanism 
of the oscillatory states initially described in Ref. [2], 
showing that the chaotic regime, being confined to far 
off-resonance cavities, is indeed mutually exclusive with 
bistable switching. 



II. MODEL DEFINITION AND LINEAR 
STABILITY ANALYSIS 
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We start from the following dimensionless coupled- 
mode model that rules the temporal evolution of the nor- 
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malized intra-cavity field a(t) coupled to the frequency 
deviation n(t), owing to the intensity-dependent refrac- 
tive index change 



da 

dt 
dn 

dt 



P + i(5 + X n ) a ~ a i 



(la) 
(lb) 



Noteworthy Eqs. describe a photonic crystal nano- 
cavity with high Q coupled to a line-defect waveguide 
P~8j [20j [24] . They implicitly assume that the nonlinear- 
ity is dominated solely by the Kerr effect with relaxation 
time r, while other possible nonlinear contributions, e.g. 
two-photon absorption along with the free-carrier disper- 
sion [HI UHl [201 E3), are neglected. Here P is the nor- 
malized power injected in the cavity through coupling 
with the waveguide, and \a\ 2 is the normalized intra- 
cavity energy, which can be easily rescaled into real- world 
units by comparison with widely used dimensional mod- 
els (see, e.g., Ref. [24]). It is worth emphasizing that 
a unit coefficient in front of the loss term in Eqs. 
implies that the time t is measured in units of the in- 
verse damping coefficient 1/ro = 2Q/wo = 2t c , where Q, 
ujq, and t c stand for the overall quality factor, the res- 
onant frequency, and the photon lifetime of the cavity, 
respectively. In these units, the two key (normalized) 
parameters are the detuning 5 = (uo — lj)/Tq, and the 
time constant r = t p Tq, where r p is the response time 
of the nonlinearity in real- world units, while % = ±1 ac- 
counts for the sign of the nonlinear Kerr coefficient. We 
point out that our model differs from [2 , inasmuch as 
the time scale is referred to the cavity lifetime instead of 
the response time of the medium. Indeed Eqs. can 
be reduced to the model analyzed in Ref. [2 by means 
of the substitution a, n, t — > a/y/r,n/r,rt. The effect of 
such transformation is however to rescale the detuning 
and the injected power in such a way that they become 
dependent on the response time of the medium, which is 
not suitable for our purpose of investigating the impact 
of the relaxation time on the dynamics of a given cavity 
with fixed characteristics. 

For a cw driving P = constant, Eqs. have the fol- 
lowing steady-state solution a(t) = A, n(t) = N = \A\ 2 , 
where 



P = E [(1 + (S + X E) 2 ] 



(2) 



E = \A\ 2 being the stationary intra-cavity energy. It 
is well known that bistability occurs for 5 > V^3 when 
X = — 1, and 5 < —V3 when \ — 1 [14 . In the discussion 
below, we will focus on the latter case, where the cavity 
resonance is blue-shifted due to the nonlinearity, a case 
which is directly comparable with the net effect of free- 
carrier dispersion induced by two-photon absorption [13] . 
All the conclusions of this paper remain valid also for x — 
— 1, provided 5 —> — 5. The values of intracavity energies 
corresponding to the knees of the bistable response are 



-2 X 5 ± V5 2 -3 



and the corresponding input powers P 6 = P(E^) can 
be calculated by means of Eq. (|2). 

The stability of the solution (2l can be investigated 
by plugging into Eqs. the ansatz a(t) — A + 5a(t), 
n(t) = N + Sn(i), while retaining linear terms in the 
perturbations 5a, Sn. 

The perturbation column array e = (5a, 5a* ,5n) T is 
found to obey the following linearized equation 



de 
dt 

M 



Me: 
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where 5 = 5 + x^- 

The characteristic equation of M reads as 

A 3 + a 2 A 2 + aiA + a = 0, 



(5) 



2 + ^a 1 = (l + ^ 2 + f), 



where the coefficients are a 2 = 

and a = ± (l + 5 2 + 2 X E5^j . 

SP occurs when the system undergoes a Hopf bifurca- 
tion, i.e. a pair of complex conjugate eigenvalues A#±iA/ 
crosses into the right half complex plane, entailing an ex- 
ponential growth of a pulsating perturbation with period 
T = 27r/|Aj|. The bifurcation point (Xr = 0) corresponds 
to the constraint a\a^ = a§, which can be solved to yield 
the following explicit expression for the SP (Hopf) thresh- 
old values E^ 
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(6) 
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and the corresponding injected power threshold Pjj = 
E±[(1 + (S + X E±) 2 }. 

The analysis reported above shows that the bistable 
response depends only on the detuning, while the time 
constant r can affect qualitatively the onset of SP due 
to the Hopf bifurcation. In fact different scenarios are 
possible depending on the existence of one or both Hopf 
thresholds (in turn corresponding to the roots in Eq. (|6|) 
being real) , and whether such threshold occurs at powers 
above or below the bistable knee for up-switching. Four 
possible scenarios are displayed in Fig. [l] where we show 
the bistable stationary response along with the unsta- 
ble eigenvalues responsible for instabilities. First, Fig. [l] 
shows the well-known fact that a purely real and pos- 
itive eigenvalue of M leads to instability of the steady 
state solution along the negative slope branch, which 
turns out to be a saddle point in phase-space. Con- 
versely, SP is characterized by a pair of complex con- 
jugate eigenvalues, and the relative threshold are high- 
lighted by empty triangles. We find such threshold to 
occur always on the upper branch of the bistable response 
(or when the response is monotone, see below). However, 
depending on the value of r, the lower Hopf threshold 
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FIG. 1. (Color online) Steady-state response E vs. P for 5 = 
—4 (a-b-c) and 6 = — 10 (d) and different r (% = 1). Stable 
and unstable branches are reported as solid and dotted lines, 
respectively. The blue and red curves superimposed on the 
right side show how the real part Re(X) and imaginary part 
Im(X) (of the dominant eigenvalue underlying the instability) 
change with E. Bifurcation points are highlighted over with 
the dashed line Re(X) = 0. The shaded regions labeled Bl 
(light blue) and SP (light yellow) correspond to the negative 
slope branch of the bistable response (real eigenvalue) and 
SP instability (pair of conjugate eigenvalues with positive real 
part), respectively. Four different scenarios are shown: (a) SP 
occurs at P = P^ values above the bistable knee value P 6 + , 
and is unbounded for increasing P; (b) SP occurs at P = P# 
values below the bistable knee P 6 + , still being unbounded; 
(c) as in (a), except SP occurs in a finite range below a given 
value P = P+; (d) as in (b), except SP occurs in a finite range 
below a given value P = P+. The shaded green regions yield 
the range of power where bistable up-switching to a stable 
state is permitted, while red ones identify the coexistence of 
a SP and an unstable saddle branch. 



P(E H ) can take place above [as in Fig. |l[a,c)] or be- 
low [see Fig. [ljb,d)] the knee P(E^) which characterize 
the bistable up-switching. In the former case the cav- 
ity can exhibit bistable up-switching to a stable steady 
state, whereas in the latter case up-switching occurs in- 
evitably towards a SP-unstable state: stabilization re- 
quires to move along the hysteresis cycle below the Hopf 
threshold. Moreover, also depending on the value of r, 
the upper branch can be indefinitely SP-unstable above 
the first threshold P(E^) [see Fig. ma, b)], or, viceversa, 
exhibit a SP switch-off energy or secondary threshold be- 
yond which SP ceases to take place [see Fig.[TJc,d); note 
that we change the detuning in Fig. [IJd) only to make 
the picture clearer, though the same qualitative behavior 
occurs at 5 = —4]. Indeed, if r > 1, Eq. (|6| may also 
admit a real solution E~^, and hence SP occurs only in 
a finite range of energies (and powers) E^ < E < E^. 
When r grows, this finite interval shrinks, and vanishes 




FIG. 2. (Color online) SP threshold energies E^ as a function 
of r for fixed detuning 5 — —4. The shaded area corresponds 
to the domain Ejj < E < E^ where SP occurs, which lies 
above the upper knee level of energy E^ (red dashed line). 




FIG. 3. (Color online) Color level plots of the (a) "on" E~ 
and (b) "off' E~^ values of threshold energy for SP, in the pa- 
rameter plane (r,S). Bist ability occurs below the line 5 = 5b. 
The curve P(E^) = P(E b ~) (black solid) divides the bistable 
region into a domain labeled BI+SP where SP sets in only 
for powers above the bistable knee for up-switching, and a 
domain labeled SP, where stable up-switching is not possible, 
being hampered by SP, which dominates the dynamics. 



as a critical value r = r c is reached. For r = r c , SP no 
longer takes place and the entire upper branch becomes 
stable. From Eq. ([6]), we find such critical value r c to be 
given by the positive real root of the cubic polynomial 
(its explicit expression is too cumbersome) 



1 + — 



1 = 0, T C 



\S\ » I- (T) 



The behavior discussed above can be clearly seen by 
reporting the SP threshold energies E^ versus r, at con- 
stant detuning. Such plot, displayed in Fig.[2]for S = —4, 
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shows a shaded region (E~ < E < £"+), which corre- 
spond to SP. We clearly see that no SP occurs for r > r c . 
Conversely, decreasing r below r c results into widening 
the portion of the upper branch that exhibits SP, until 
below r = 1 the whole upper branch above the switch- 
on threshold E]j turns out to be unstable. In this case, 
the SP switch-off energy E^ diverges as the asymptote 
r = 1 (dashed vertical line) is approached. Importantly 
for r — >• also the first threshold E~^ diverges, which 
means that a finite response time is a key ingredient for 
SP to be observable. This is consistent with the fact that 
Kerr instantaneous nonlinearities yield no SP at all. In 
fact, in this c ase, the eigenvalu es are easily found to be 
A ± = -1 ± y/E 2 - (5 + 2xE) 2 , which rule out the possi- 
bility to have a complex conjugate pair with positive real 
part. 

In order to have a complete picture and further show 
how the onset of SP changes with detuning, we have 
drawn in Fig. [3] a color map of the level curves of SP 
threshold energies E^, in the parameter plane (r, 5). A 
number of interesting observations can be drawn. The 
SP region is bounded by the border r = r c , and in the 
bistable region (5 < 5b) r c decreases with decreasing val- 
ues of absolute detuning \5\. Interestingly enough, the 
scenario illustrated in Fig. [2] remains valid also for detun- 
ings 5 > £5, where bistability disappears. Finally in the 
region of positive detunings, we are left with the upper 
branch being fully unstable for all energies E > E~^. In 
this region, the threshold energy E~^ diverges, not only in 
the instantaneous limit r = 0, but also for r = 1, whereas 
relatively low values of E~^ are found for r ~ 1/2, i.e. 
when the response time of the nonlinearity is nearly equal 
to the photon lifetime. Furthermore, the bistable region 
is divided into two distinct domains by the (solid black) 
curve which arises from the condition P(E^) = P(E^) 
(its explicit expression is too cumbersome). In the do- 
main BI+SP above such curve (bounded from above also 
by the line 5 = 5b), one has that the cavity can work as 
a bistable switch since the upper branch right above the 
knee for up-switching is stable [as in Fig.[TJa,c)], whereas 
in the domain labeled SP below the curve, up-switching 
is no longer allowed, since the upper branch above the 
knee is SP-unstable [as shown in Fig. [TJb,d)] The reader 
can easily recognize a qualitative similarity of the picture 
discussed here with the dynamics of SP ruled by free car- 
rier dispersion, recently discussed in Ref. [13]. Although 
a detailed analytic investigation of the stability of the 
SP-oscillating state (limit cycle) is beyond the scope of 
this paper, similarly to the case discussed in Ref. [13], 
our numerical simulations of Eqs. ([!]) suggest that sta- 
ble limit cycles, working as attractors from a large basin, 
exist in a wide domain of the parameter plane (witness- 
ing the supercritical nature of the Hopf bifurcation). An 
example of such stable dynamics is shown in Fig. [4] 

Let us comment on the observability of the SP dynam- 
ics. In nanocavities with high Q (Q = 10 3 — 10 5 ) the 
photon lifetime in the near infrared ranges from few pi- 
coseconds to tens of picoseconds. While the constraint 



to have a response time of the nonlinearity of the same 
order of magnitude is naturally met in semiconductors 
with nonlinear response dominated by free-carrier disper- 
sion, the same constraint in the framework of the Kerr 
model rule out the possibility to observe SP dynamics in 
media with nonlinearities of electronic origin since they 
are too fast (fs range). Nevertheless the predictions of 
our Kerr model become interesting for Kerr-like mate- 
rials with response time in the ps range such as, e.g., 
soft matter, metal films [26], or more traditional liquids 
with reorientational nonlinearity, which are still the ob- 
ject of recent studies [27j[28]. In particular, for instance, 
highly nonlinear liquids such as nitrobenzene [29] or CS2 
could be easily employed to fill a photonic crystal ma- 
trix (as also recently proposed for microstructured fibers 
[28]). while metal films could be employed in conjunction 
with dielectrics to form a single cavity or cavity arrays 
[30] . Assuming, for instance, a response time r p h ys ~ 30 
ps [29], which yield r = (T p hy S /2t c ) ~ 0.75 in a cavity 
with Q ~ 25000 (t c ~ 20 ps at A = 1.55/im), assum- 
ing U2i ~ 10~ 17 m 2 jW and a nonlinear modal volume 
V = 3(A/n) 3 , the threshold power P = 10 in Fig. [I] 
corresponds to a real- world power Pi n = (j/Tq)P ~ 10 
mW in the waveguide coupled to the nanocavity, where 
7 = uori2ic/(n e ffnV) is the overall nonlinear coeffi- 
cient [24]. Here we have assumed a refractive index 
n ~ n eff ~ 1-5 and Q to be essentially determined by 
the coupling itself. 

Having characterized so far the threshold for SP and 
its competition with bistability, since Ikeda and Akimoto 
have shown that the limit cycles destabilize, leading even- 
tually to chaos [2], in the next section we deepen this 
point with the aim of determining the domain of the pa- 
rameter plane where the transition to chaos could be ob- 
served. 



III. THE TURBULENT REGIME 

In Ref. [2 Ikeda and Akimoto have studied the transi- 
tion to chaos, identifying a period doubling cascade up to 
2 2 P (i.e. oscillation with period-four) at a fixed value of 
detuning. Here we report further details about the emer- 
gence of chaos in a wide domain of parameters. We em- 
ploy different tools, ranging from Poicare section and its 
corresponding bifurcation diagram to the calculation of 
Lyapunov exponents. Our principal purpose is to assess 
whether a nano-cavity described by the model can 
work as a reliable bistable switch, and hence whether 
the onset of chaos should be expected when the cavity 
operates progressively off-resonance, especially in the re- 
gion labeled BI+SP in Fig. [3] To begin with, it is in- 
structive to report about the dynamics ruled by Eqs. ([!]) 
when, starting above the threshold P^, the input power 
P is adiabatically decreased. In fact, this is the situation 
where the onset of chaos is expected according to Ref. 
[2]. We start at moderately low detuning (5 = —4) and 
for r = 0.45, which corresponds to SP being unbounded 
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FIG. 4. (Color online) Dynamics of SP ruled by Eqs. (1 
lb), with 6 = —4, and r = 0.75: (a) temporal evolution 
of the intra-cavity energy and carrier density corresponding 
to the rightmost blue circle in (b); (b) steady response with 
superimposed peak energy of the periodic oscillations (blue 
open circles) ; The red filled circle marks the Hopf bifurcation 
point. The black filled circle marks the SP-unstable steady 
state which gives rise to the dynamics shown in (a,c); (c) 
phase-space picture of the optical field showing the attracting 
limit cycle from two different initial conditions (open circles) ; 
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FIG. 5. (Color online) Dynamical evolution ruled by Eqs. |T]) 
as the forcing term y/P varies adiabatically in time (dashed 
line). Here r = 0.45, and (a) S = -4; (b) 5 = -10. The ver- 
tical lines mark the time instants at which the input power 
cross the main bifurcation points: bistable knees (dotted 
blue and black, respectively) and Hopf threshold P^ (dashed 
red). The insets zoom over characteristic time intervals, in- 
dicated by arrows. Notice that in (b) the two thresholds P^ 



and P H almost overlap. 



on the upper branch. As shown in Fig. [5^a) , in this 
case, the Hopf bifurcation is clearly supercritical, since 
the system settles on a limit cycle, whose amplitude van- 
ishes approaching the bifurcation point (approximately 
as [E — Eft) 1 / 2 )' However, at larger (in modulus) detun- 
ings, we observe a jump in \a\ 2 as the limit cycle (IP) 
looses its stability and the system settles on a period- 
two (2P) oscillation, as shown in Fig. j5^b) for 8 — — 10. 
The 2P solution does not visit anymore the simplest IP 
limit cycle. Conversely it abruptly switches to the sta- 
ble low-branch steady state. Remarkably this happens 
still above the Hopf bifurcation point E~^. A more com- 




FIG. 6. (Color online) As in Fig. [5] with S = -15. The 
chaotic region is highlighted in yellow. 



plex sequence of period doubling bifurcations and chaotic 
motion are detected at higher detunings when the input 
power approaches the knee value P 6 ~ . In Fig. [6j where 
5 = —15, oscillations with several different periods are 
evident. Moreover a chaotic regime appears, in the range 
\[P « 17 - 20 (P « 280 - 400). In phase space this cor- 
responds to the appearance of a strange attractor (not 
shown because its structure is already illustrated in Ref. 

0). 

Since following simply the adiabatic dynamics could 
be possibly misleading (e.g. because of critical slowing 
down), as it neglects the rich variety of phenomena that 
occurs over the small scale, we have drawn also bifurca- 
tion diagrams calculated by collecting trajectory points 
on a Poincare section for different powers. A typical ex- 
ample, using the same parameters as in Fig.[6| and defin- 
ing the Poincare section on the fixed phase Za — tt = 
of the intra-cavity field, is reported in Fig. [7| We can 
clearly identify a period doubling cascade (up to 2 3 P), 
as well as chaotic regimes. The onset of chaos follows 
a non-trivial scenario where narrow windows of period- 
three solutions (3P) are interspersed between two ranges 
of powers where the motion turns out to be chaotic. The 
vertical dashed lines in Fig. [7] mark indeed a 3P window. 
This is analogous to the Feigenbaum's route to chaos and 
confirms the observation of chaos for P « 220 — 380, al- 
ready drawn above from Fig. [6] 

The bifurcation diagram is computed up to P « 150 
because lower power levels do not result into any limit 
cycle. Vice-versa the solution is observed to collapse to- 
ward the stable node represented by the lower branch 
of the bistable response. This phenomenon is indepen- 
dent from chaos, as mentioned above with reference to 
detuning S = —10. It occurs at large negative detun- 
ings for input powers above the SP threshold (P > P#). 
This can be qualitatively explained by the coexistence 
of a stable fixed point (lower branch solution), a saddle 
(negative-slope branch), and an unstable limit cycle. As 



FIG. 7. (Color online) Bifurcation diagram, r = 0.45, 6 = 
— 15. The vertical dashed lines highlight a 3P window and 
the collapse at small P. Inset shows a detail in which period 
doubling and periodic windows can be identified more clearly. 



the 2P limit cycle spans the phase space with wide os- 
cillations around the upper branch, it can approach the 
saddle point (which near the first bistable knee is closer 
in phase space to the center of the oscillations) being 
forced away, until eventually it can be captured by the 
lowest energy (stable) solution. 

While the bifurcation map is a useful visual tool to 
characterize the onset to chaos and the full dynamics at 
fixed parameters, in order to explore in which region of 
the parameters one should expect to observe the chaotic 
dynamics, we have resorted to compute the dominant 
(maximal) Lyapunov exponent. We have explored a wide 
region of the parameter plane (r, 5), where, in each point 
of such plane, we have iterated over the values of power 
P to find the largest exponent. We recall that a posi- 
tive Lyapunov exponent (within numerical inaccuracies) 
entails that the system exhibits a chaotic behavior (here 
quasi-periodic motion is excluded by the dissipative char- 
acter of our model (la)). From the map displayed in 
Fig. [8j we notice that chaotic motion manifests itself only 
when the cavity is detuned far off-resonance (i.e. at very 
large values of \S\), provided that the SP-unstable range 
is not finite, or in other words that the Hopf bifurcation 
is not bounded from above (E^ — >• oo) which requires 
r < 1. Therefore we can draw the important conclu- 
sion that the region where the cavity could work as a 
bistable switch is mutually exclusive with chaos. There- 
fore the onset of chaos cannot spoil the behavior of the 
cavity as a switch, once the latter is used in the region 
labeled BI+SP in Fig. [3j We point out that, in terms of 
power, the observation of chaos is much more challenging 
than stable SP since power levels leading to the former 
turn out to be much larger than those leading to the lat- 
ter; indeed at very large detuning, which corresponds to 
several times the cavity linewidth, bistability is observed 
at much higher power level (compare the horizontal axis 



FIG. 8. (Color online) Color level map of maximal Lyapunov 
exponent in the parameter plane (r,S). 



scale in Fig. [4] and Fig. [7]). 

As a final remark about the existence of 2 n P periodic 
solutions, we point out that they can be detected only 
when Pft < P 6 ~ [as in the examples shown in Fig.[TJb,d)]. 
As discussed with reference to Figs. [T]|3j this may occur 
not only for both r < 1 (as shown explicitly above), but 
also for r > 1, where the Hopf bifurcation is bounded 
from above. More specifically 2P solutions are easily seen 
in a small subset of the region marked as SP in Fig. [3] 
In this case two unstable solutions, namely a repulsive 
(negative slope) branch and a SP branch, coexist. This 
seems to be a key ingredient for the limit cycles to loose 
their stability. 



IV. CONCLUSIONS 

In this work we have revisited the model that rules 
the behavior a passive small cavity with Kerr delayed re- 
sponse, pioneered in Ref. [2 . We have reported a full 
characterization of SP instabilities and their competition 
with bistability, outlining the existence of different pos- 
sible scenarios. Importantly we have found a maximal 
critical value for the relaxation time that allows SP to 
occur, and have shown that SP can have two bifurcation 
points, while it can occur also in the absence of bista- 
bility. We have further characterized the destabilization 
mechanism of the limit cycle in the full parameter space, 
finding that chaos is mutually exclusive with the domain 
where the cavity can be employed as a bistable switching 
element. In particular the chaotic regime predicted by 
Ikeda [2 in this system requires to be detuned strongly 
off-resonance, in turn implying the use of extremely high 
powers, thus making its observation rather challenging. 
Viceversa, in contrast with Kerr instantaneous nonlin- 
earities, the observation of stable SP appear feasible in 
high-Q nanocavities filled with Kerr-like media with re- 
sponse time in the range of picoseconds. Future work 
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will be devoted to study the effect of coupled cavity sys- (no. 249012). 
terns and the interplay of different nonlinear mechanisms. 
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